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Abstract 

Very narrow spatial bright solitons in (1+1)D and (2+l)D versions of cubic-quintic 
and full saturable models are studied, starting from the full system of the Maxwell's 
equations, rather than from the paraxial (NLS) approximation. For the solitons with 
both TE and TM polarizations, it is shown that there always exists a finite minimum 
width, and they cease to exist at a critical value of the propagation constant, at which 
their width diverges. Full similarity of the results obtained for both nonlinearities 
suggests that the same general conclusions apply to narrow solitons in any non-Kerr 
model. 



1. Introduction 

The standard approach to description of spatial solitons in nonlinear optical waveguides is 
based on the use of the nonlinear Schrodinger (NLS) equation, which replaces the Maxwell's 
equations (ME) in the paraxial approximationEI (see also Refs.@). As is known, this approxi- 
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mation can be insufficient for the description of very narrow, subwavelength (subwavelength) 
solitons, with a size ^ A, where A is the carrier wavelength in vacuumi. Terms which are 
neglected when deriving the NLS equation from ME couple to the propagation constant, 
affecting dynamics of very narrow solitons with arbitrary polarizationi. Particular cases of 
the TM (transverse-magnetic)- and TE (transverse-electric) polarized subwavelength soli- 
tons in media with purely cubic (Kerr) nonlinearity were considered in detail in Refs.ii. 
It was found that, in the TM case, the size of both bright and dark solitons cannot be 
essentially smaller than A/2. In the TE case, the size of the dark soliton is also limited 
from below, while the bright soliton may formally be arbitrarily narrow; however, narrow 
TE solitons are subject to a strong instability. 

The analysis performed in those works was formal in the sense that subwavelength soliton 
solutions in the model with the Kerr nonlinearity imply unrealistically large values of the 
nonlinear correction An to the refractive index. An ^ 1. As is knownS, the quadratic (in the 
amplitude of the electromagnetic field) correction An, which defines the Kerr nonlinearity, 
IS, as a matter of fact, only the first term of the expansion in powers of the squared field. 
A detailed analysis shows that results of the formal consideration, based on dropping the 
higher-order corrections and using the ensuing truncated model in the range where the 
quadratic correction An is allowed to be large, are ambiguous: they strongly depend on a 
particular stage of the analysis at which the higher-order terms are omitted. This ambiguity 
is especially conspicuous in the case of the TM narrow solitons: depending on the choice of 
the truncation stage, one arrives at a conclusion that the width of the bright TM soliton 
remains limited from below, or may become arbitrarily small. 

In fact, the physical problem of the existence of subwavelength solitons is unresolvable 
within the framework of the Kerr model, and the only possibility to obtain a definite result 
is to adopt a more realistic nonlinearity, with saturation of An, in some form for strong 
fields. While an exact form of the saturation in real materials is not usually known, there 
are two commonly used models, viz., the saturable nonlinearity properi, and the cubic- 
quintic (CQ)0i one. The former nonlinearity is realistic for alkali metal vaporsB, and the 



latter nonlinearity, which combines self- focusing cubic and seli- defocusing quintic terms, is 
well documented to describe nonlinear optical properties of the PTS crystalS. Moreover, 
formation of (2+l)D solitons in this material, consistent with the predictions of the CQ 
model, was recently observed in the experiment reported in Ref.0. We also note that the 
(2+l)D and (3+l)D versions of the CQ model have recently attracted considerable attention 
as they give rise to stable or almost stable multidimensional (spatial or spatiotemporal) 
solitonJi. 

Our objective is to consider subwavelength spatial solitons of both TM and TE types 
parallel in the CQ and saturable models, in order to arrive at general conclusions con- 
cerning their width and other fundamental properties. Moreover, both (1+1)- and (2+1)- 
dimensional [(1+1)D and (2+l)D] cases will be considered, in the latter case the spatial 
soliton being a cylindrical beam. We will conclude that, in both models and for both dimen- 
sions, there always exists a finite minimum value of the soliton's width, and fundamental 
properties of the subwavelength solitons are fairly similar in these two models. On the basis 
of these results, one may infer that, in fact, qualitative results obtained in this work (first of 
all, the existence of the minimum width) apply to narrow spatial solitons in any non-Kerr 
model. 

The rest of the paper is organized as follows. In section 2, we derive, directly from ME, 
general equations that describe narrow spatial TM solitons in both (1+1)D and (2+l)D 
geometries. In section 3, we consider a much simpler TE case. In fact, for this case soliton 
solutions are well known (in particular, they can be found in an exact form in the (1+1)D 
geometry). This allows us to rigorously prove that the relative soliton's width, divided 
by the carrier wavelength, cannot be smaller than a certain minimum value. In section 
4, combining analytical considerations and direct numerical solutions, we find (1+1)D and 
(2+l)D fundamental TM solitons and arrive at the conclusion that their width is also limited 
from below. 
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2. Equations for the transverse-magnetic case 

We start the analysis with the most nontrivial case of the TM polarization. In this case, 
irrespective of the particular nonlinearity, the real (physical) vectorial electric £ and magnetic 
Ti fields can be taken as 

TLy = -H{x) cos n-^ = n^ = 0, (1) 

So: = E{x) cos$, £z = E^{x) sin$, £y = (2) 
in the (1+1)D case, and in the (2+l)D case, 

= {y/r)H{r)cos<^, Hy = -{x /r)H{r) cos^, = 0, (3) 



£x = (x/r)£'(r) cos$, £y = {y /r)E{r) cos ^, E-^ = i?^(r) sin $; (4) 

note a phase shift 7r/2 between the transverse and longitudinal components in these expres- 
sions. Here, x,y and z are the transverse and propagation coordinates, = x"^ + y^, t is 
time, and the common phase of all the fields is $ = jSz — ujt, (3 and uj being, respectively, the 
propagation constant and frequency of the carrier wave. In the (2+l)D case, the expressions 
(0) and (0) correspond to the standard TM^i mode in a cylindrical waveg uiddi. 

The simplest saturable model assumes an isotropic material characterized by the relation 



1 + 



(5) 



between the electric induction and strength, with positive nonlinear permeabilities £2 and 
£4. Expansion and truncation of Eq. (|^) yields the CQ model in the form 

PcQ = ^ (^0 + ^2^' - e^£^) . (6) 

A difference between the two models appears in the case when the truncation leading 
from Dgat to T>cQ is no longer valid. Below, we present analysis in a detailed form for the 
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CQ model. For the saturable one, it is quite similar, but formulas are more cumbersome. 
Final results will be displayed for both models together. 

Following the usual rotating-wave approximation, the next step is to substitute Eqs. (|I]) 
and (0) or and (H) into the relation (|^) (or (|^), in the case of the saturable nonlinearity) , 
and collect all the contributions to the fundamental harmonics sin<l> and cos$, neglecting 
higher-order harmonics. This yields direct relations between the electric-field induction and 
strength, 

where, in the case of the CQ nonlinearity, effective nonlinear susceptibilities are 

et^e, + \e2{^E^ + El) - ^6,i5E^ + 2E'eI + E^), (8) 

ei = eo + -ME'' + ?,El) - \e^{E^ + 2E''EI + bEt). (9) 
4 8 

Insertion of the above expressions and @ or (^) for the magnetic and electric field 
and (^) into the Maxwell's vectorial equation for the electric field, 

^ = 

C Ot 

where c is the light velocity in vacuum, we arrive at a scalar ODE 

E[+PE = -{uj/c)H, (10) 

where the prime stands for d/dx or d/dr in the (1+1)0 and (2+l)D cases, respectively. 
Further, we substitute the (1 + 1)D or (2+l)D expressions ([^) and (Q) or (^) and (^), in 
combination with the relations (|^) and (H), (^) which define the electric induction, into the 
Maxwell's vectorial equation for the magnetic field, 

c ot 

This yields two more equations, one of which is again an ODE, while the other one is just 
an algebraic relation. 



^^^i-D(^D-i^y = {u:/c)E,eu (11) 
(3H=-{u/c)E, (12) 

where Z) = 1, 2 is the transverse dimension. 

Ehminating the magnetic field H from Eqs. ([10|) and (|ll|) by means of Eq. (p!2D, we then 
obtain a system of two equations for the fields E and Ez, 

K = E{l-P), (13) 
^l-D (rD-iEvi)' = -EzV2. (14) 

Here, the variables have been rescaled as {x,y,r) — > Py/^{x,y,r), [E,Ez] —>■ 
E, {EzP/y/^)\ {uj/p^)y/3£^/2, with 7 = 1- {Po/Pf, where po = y/e^uj/cis the propaga- 
tion constant in the linear regime, and 

P^E^ + IeI + a[E^+ ^-Ie^EI + :^Et^ , (15) 
^1 = 1 + (16) 



V2 = l + 
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Ie^ + ^El + a i^-E' + ^E^El + ^^E^ 



(17) 



1-7 

a = -7(l-7)''s, (18) 
s = {10/9)6o£4e2^. (19) 

These equations will be used below to study TE and TM solitons in the CQ model, 
which will be paralleled by the same analysis for the model with the saturable nonlinearity. 
Before proceeding to that, it may be relevant to revisit the case of the TM soliton in the 
Kerr (rather than CQ) (1+1)D medium. In that case, it is necessary to bear in mind the 
pure cubic (Kerr) nonlinearity applies as long as the nonlinear correction to the refractive 
index is much smaller than the linear index. In the present notation, this condition can be 
shown to amount to inequalities 

E\x) < land {Ef < 1. (20) 
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Taking these into regard and performing the corresponding expansions in the above 
equations (p^ through (p]7|), one can ehminate the longitudinal electric field and derive 
an eventual equation for the transverse field, 



E" --fE + 



E^ + - {E 



/\2 



E = 0, 



(21) 



where 7 was defined by Eq. (??)• In comparison with the traditional paraxial (NLS) 
approximation, the only new term in Eq. (|21|) is (1/3) {E')'^ E, which, as a matter of fact, 
is a contribution from the longitudinal component E^ in the present TM case. 



Obviously, Eq. (pT]) can give rise to bright solitons only in the case 7 > 0, which will 
be assumed to hold hereafter. Note that the classical broad solitons correspond to the case 
7 <^ 1. In that case, the broad soliton with the first correction produced by the extra term 
can be easily found: 



Esoi(x) 



1 
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(22) 



^l + _Ljsech(V70-^sech^(V70 • 

Numerical solution of Eq. (^) shows that it does give rise to solitary-wave solutions 
with an arbitrary small width; however, when the soliton becomes too narrow, the solution 
violates the applicability conditions (pOD, which makes it necessary to modify the nonlin- 
earity, i.e., to consider the model with the CQ or saturable nonlinearity, which will be done 
below. 



3. The transverse-electric case 

The above equations, derived for the CQ nonlinearity, also contain an essentially simpler 
case of the TE polarization, which can be obtained setting formally 7 = but keeping a 7^ 0. 
In this case, vi = V2 = 1-, P = E"^ + aE^, and in the (1+1)D geometry an eventual equation 
for E{x) takes a well-known form 

E" = E- E^ - aE^, (23) 

a commonly known exact soliton solution to which can be written as 
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4/E\x) = 1 + 1^1 + 16a/3cosh(2x). (24) 
It exists provided that a > —3/16, or, with regard to Eq. (|T^, 

(/?//?o)' < (m)cr = l + 3/(16s). (25) 



As {P/PoY is approaching the value {(3 / [3q)1^, the average width of the sohton (^41), defined as 
per Eq. (|30D, diverges as |ln [1 + 3/(16s) - {l3/(3of]l while the soliton's amplitude remains 
finite, Emax = 2. 

The main characteristic of the spatial soliton is its width. For the (1+1)D soliton (p4D, 
the FWHM0 width A can be easily found, 

cosh A = (1 + 16a/3)"^/^ + 2. 

In unnormalized units, the FWHM width is 

W = (4A/37rno)VsF(16a/3), 

where 

F{y) = l^r^^^ln 



1 + y)~'^^ + 2 + J3 + 4(1 + y)-i/2 + (1 + y)-i 



(26) 



The function F, and hence the width of the TE soliton, attain a minimum value at cr = 
—0.1493. The ratio of the corresponding minimum width of the TE soliton, normalized to 
the carrier wavelength, is 

wZn^/X = 0.825v/^/£2. (27) 

In the next section, we will also use a different (average) definition of the width, based 
on Eq. (pOl). It is easy to find that the minimum value of the average width (|30|) differs 



from the minimum FWHM value (pT]) by less than 8%. 

Consideration of TE solitons in the (2+1 )D geometry leads to an equation 

d'^E IdE 1 ^ ^ . . 

-— + -E = E - E^ - aE'' 28 

ct^r r ar 



(r again being the radial variable), cf. Eq. (pS]). In fact, Eq. (^) is exactly the same 
equation which describes known solitons with an internal spin (vorticity) s = 1 in the CQ 
model in the (2+l)D geometryH. The origin of an effective "vorticity" term, — i?/r^, in Eq. 
(p8|) is the structure of the (2+l)D ansdtze @ and (|^) for the vectorial fields. 

Equation ( pS] ) has been studied in detail numerically elsewhereil. Although the minimum 
width of a family of solitons generated by this equation was not specially considered, a 
straightforward consequence of results presented in Refs.i is that the minimum width is 
always finite. 

Thus, it is possible to rigorously prove the existence of a finite lower bound for the width 
of spatial TE solitons in the (1+1)D and (2+l)D models with the CQ nonlinearity. Quite 
similarly, the same results can be obtained for TE solitons in the model with the saturable 
nonlinearity (we do not display technical details here, as they do not contain anything 
essentially novel). 



4. Transverse-magnetic narrow spatial solitons 

The most interesting case is that with the TM polarization, as it does not reduce to a 
single-component equation. As it was mentioned above, the necessary condition for the 
existence of localized solutions to the corresponding equations (|l^) and (p^ , supplemented 



by Eqs. (|T^) through ([T9|) , is (3q < (3"^. Below, we will use a relative propagation constant, 
P/Pq = (1 — 7)~^/^, as a measure of the departure from the paraxial approximation: in the 
limit of a very broad soliton, one has f3/(3o 1; while in the opposite limit of an infinitely 
narrow soliton (if any), P/Po — oo. 

A straightforward analysis of Eqs. ([T3|) and ([M]) makes it possible to prove that, in the 
(1+1)D case for both CQ and saturable nonlinearities, localized solutions with no zeros of 
E{x) and Ez{x) at a; 7^ do not exist if £'(0) = 0. Further analysis, details of which are not 
displayed here, shows that, in the (1+1)D geometry, a localized fundament al-soliton solution 
has an even monotonically decreasing transverse component ii^(|x|), while the longitudinal 
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one Ez{x) is odd, with Ez{0) = and a single extremum at finite see Fig. la (in this 
and next figures, the results are simultaneously shown for the CQ and saturable models). 
In the (2+l)D case, fundamental solitons feature a different structure: they have E{0) = 0, 
and both E{r) and Ez{r) may have zeros at finite r, see Fig. lb. 

Proceeding to the width of the TM soliton, we use the electromagnetic energy density. 



The FWHM definition of W, which was used in the previous section for the TE soliton, 
is ambiguous for TM solitons, as they have three different components. The results for 
the average width, which were obtained from numerical solutions of Eqs. (pTS]) and (^4|), 
are summarized in Fig. 2, showing W vs. the relative propagation constant at different 
values of the material constant s defined by Eq. (^). Basic features of the dependence 
are the same as those following from the above analytical expressions for the TE soliton 
in the (H-l)D case: there is a finite minimum of the soliton's width, and the solitons do 
not exist beyond a critical value of f3/(3o, at which the soliton's width diverges while its 
amplitude remains finite, cf. Eq. (|25|) . Naturally, the minimum value of the width depends 
on the dimension of space, type of the nonlinearity, and, for a fixed nonlinearity, on values 
of material constants. As well as in the case of the TE solitons, the minimum width tends 
to zero and the maximum value of P/Po diverges as £4 — > 0. 

Note that the soliton shapes shown in Fig. 1 pertain to the fixed value of the relative 
propagation constant f3/(3o = 1-05, which is rather close to the paraxial limit. Comparison 
with other numerical results demonstrates that, for both CQ and saturable models, the 
shapes of the (H-l)D and (2+l)D TM solitons remain fairly similar to those shown in Fig. 
1 with the increase oi (3/ (3q up to the point where the minimum width is attained. With the 
subsequent increase of /?//?o up to the value at which the solitons cease to exist, the shape 
changes much more. 



(29) 



to define the soliton's average half-width as 




(30) 
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The above considerations did not tackle the stabihty problem. While a consistent sta- 
bility analysis requires very tedious simulations of the Maxwell's equations with the CQ or 
saturable nonlinearity, which is beyond the scope of this work, we tried to test the stability 
by means of the simple Vakhitov-Kolokolov (VK) criterionS (although its applicability to 
the present model is not obvious): the propagation constant must have a positive slope as 
a function of the beam's power (norm of the solution), 



(cf. Eq. (0)), which is a dynamical invariant of the model. Our numerical calculations 

have shown that all the TM solitons, in both the (1+1)D and (2+l)D geometries, satisfy 
the VK criterion. 

5. Conclusion 

In this work, we have studied very narrow (1 + 1)D and (2+l)D spatial TE and TM soli- 
tons supported by the quintic-cubic and saturable nonlinearities, starting the full Maxwell's 
equations. For TE solitons, we have obtained a single equation, which is equivalent to the 
stationary version of the nonlinear Schrodinger equation with the corresponding nonlinear- 
ity, while for TM solitons we end up with a system of equations for the transverse and 
longitudinal components of the electric field, rather than with a single one. It is also note- 
worthy that TE solitons in the (2+l)D case are described by the same equation as vortex 
solitons with spin s = 1 in the (2+1 )D nonlinear Schrodinger equation. 

Detailed analysis was displayed for the cubic-quintic model, while final results were 
given for the saturable model as well. Full qualitative similarity of the results obtained for 
both models (although exact scales may be quite different, see caption to Fig. 2) makes it 
very plausible that the situation is basically the same for any physically relevant non-Kerr 
medium. The most important conclusions are that, in all the cases, there is a finite minimum 
of the soliton's width, and the solitons cease to exist at a critical value of the propagation 
constant, at which their width diverges. 
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Figure captions 



Fig. 1. Typical examples of the fundamental spatial TM soliton in the (1+1)D (a) and 
(2+l)D (b) geometries. The relative propagation constant is /5//?o = 1-05, and the material 
constant (|l^) is s = 1. In this and next figures, the solid and dotted curves pertain, 
respectively, to the cubic-quintic and saturable models. 

Fig. 2. The width of the (1+1)D (a) and (2+l)D (b) fundamental spatial TM solitons 
vs. the relative propagation constant (3 / (3o at various fixed values of the material constant 
s. For the saturable model, the deviation of the relative propagation constant from 1 is five 
times that shown on the horizontal axis for the cubic-quintic model. 
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